Spin 3/2 Penta-quarks in anisotropic lattice QCD 
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A high-precision mass measurement for the pentaquark (5Q) O 1 



J 1 



3/2 ± channel is 



performed in anisotropic quenched lattice QCD using a large number of gauge configurations as 
N con f = 1000. We employ the standard Wilson gauge action at [3 = 5.75 and the O(a) improved 
Wilson (clover) quark action with n = 0.1210(0.0010)0.1240 on a 12 3 x 96 lattice with the renor- 
malized anisotropy as a s /a t = 4. The Rarita-Schwinger formalism is adopted for the interpolating 
fields. Several types of the interpolating fields with isospin 1 = are examined such as (a) the 
NK*-type, (b) the (color-)twisted NK*-type, (c) a diquark-type. The chiral extrapolation leads to 
only massive states, i.e., msQ ~ 2.1 — 2.2 GeV in J p = 3/2~ channel, and msq = 2.4 — 2.6 GeV in 
J F = 3/2 + channel. The analysis with the hybrid boundary condition(HBC) is performed to inves- 
tigate whether these states are compact 5Q resonances or not. No low-lying compact 5Q resonance 
states are found below 2.1GeV. 

PACS numbers: 12.38.Gc, 12.39.Mk, 14.20.-c, 14.20.Jn 
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I. INTRODUCTION 

The recent discovery of the manifestly exotic baryon 
6+ (1540) by the LEPS group at SPring-8 has made a 
great impact on the exotic hadron physics pj. Apart 
from other pentaquark baryon candidates, S~~(1862) 
2] and 9 C (3099) @, several other candidates of exotic 
hadrons have also been discovered, such as A" (3872), 
A,(2317), 50(3115), AT(3940) and F(3840) A\. They 
also receive an increasing interest from theoretical side 
as well. G + (1540) is supposed to have baryon number 
B = 1, charge Q = +1 and strangeness S = +1. Since 
uudds is the simplest quark content to implement this 
quantum number, + (154O) is a manifestly exotic penta- 
quark(5Q) state. The penta-quark + had been consid- 
ered several times even before the experimental discov- 
ery H H HH 0. In particular, Ref.@ provided the 
direct motivation of the experimental search [|. The 
discovered peak in the nK + invariant mass is centered 
at 1.54 ± 0.01 GeV with a width smaller than 25 MeV. 
At the present stage, some groups confirmed the LEPS 
discov ervllOl IrH Il2l Il3| , while the others reported null 
results 14]. It will still take a while to establish the exis- 
tence or non-existence of O + (1540) experimentally [l5| . 
Ref.^3 claims that + must be isoscalar, since no ++ 
is observed in the pK + invariant mass spectrum. 

Enormous theoretical efforts have been devoted to 50 
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One of the most challenging problems in understanding 
its structure is its extremely narrow decay width as T < 1 
MeV [55|. Several ideas have been proposed: (1) I = 2 
possibility |§l|, (2) Jaffe-Wilczek's diquark picture [22| . 
(3) ttKN hepta-quark picture [2E IH3. (4 ) string picture 
go, EH, (5) J p = 3/2" possibility^ (H. Although 
each gives a mechanism to explain the narrow decay 
width, none of them can satisfy all the known proper- 
ties of O + (1540) simultaneously. 

In this paper, we are interested in the J = 3/2 pos- 
sibility, J p — 3/2~ in particular. Note that the spin 
of O + (1540) has not yet been determined experimen- 
tally. In the constituent quark picture, the narrow decay 
width of J p = 3/2~ penta-quarks can be understood 
in the following way 0, E^] . We expect that the spe- 
cial configuration (Os) 5 is dominant in the 5Q ground- 
state in J p = 3/2~ channel. Although J p = 3/2 _ 
penta-quarks can decay to KN in the d-wave, the spec- 
troscopic factor to find d-wave KN states in the dom- 
inant (Os) 5 configuration vanishes. Since the decay is 
thus allowed only through its sub-dominant d-wave con- 
figuration, the decay width is suppressed. Note that 
it is further suppressed by the d-wave centrifugal bar- 
rier, leading to the significantly narrow decay width of 
J p = 3/2~ penta-quarks. A possible disadvantage of 
J p = 3/2~ assignment is that such a state tends to 
be massive due to the color-magnetic interaction in the 
constituent quark models, which seems to be one of 
the main reasons why there are only a limited num- 
ber of effective model studies for spin 3/2 penta-quarks 
[HElElElEEElElEiEElHlElEi However, 
it is not clear whether these conventional framework is 
applicable to a new exotic 5Q system as O + (1540) with- 
out involving any modifications. Indeed, a model was 
proposed where a part of the role of the color-magnetic 
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interaction can be played by the flavor-spin interaction, 
which makes the mass-splitting between the 1/2 - and 
the 3/2~ states smaller (45j. 

There have been several lattice QCD calculations of 5Q 
states by today [H IH IH HJ |H M |H |H HI il |^ . 
However, these studies are restricted to J p — l/2 ± chan- 
nels except for a very recent one |67|. Enormous efforts 
are being devoted to more accurate studies of J p = l/2 ± 
states, using the variational technique to extract mul- 
tiple excited states, among which a compact resonance 
state is sought for. Indeed, quite larg e scale calculations 
are planned and being performed 68] attempting to elu- 
cidate some of the mysterious natures of O + (1540) such 
as its diquark structure and/or non- localities desired in 
interpolating fields. Here, we emphasize again that these 
studies are aiming at J p = l/2 ± states, not at 3/2 ± 
states. 

In this paper, we present anisotropic lattice QCD re- 
sults on 5Q states in J p = 3/2 ± channels using a large 
number of gauge configurations as N con f = 1000 as an 
attempt to search for a low- lying 5Q state in J p = 3/2 ± 
channel. We adopt the standard Wilson gauge action 
at j3 = 5.75 on the 12 3 x 96 lattice with the renor- 
malized anisotropy a s /a t — 4. The anisotropic lattice 
is known to serve as a powerful tool for h igh- precision 
measurements of temporal correlators [6|| [qJ, El E3- 
The large number of gauge configurations -/V con f = 1000 
plays a key role in our calculation, because 5Q corre- 
lators in J p = 3/2 ± channels are found to be quite 
noisy. For quark action, we adopt 0(a)-improved Wilson 
(clover) action with four values of the hopping parameter 
as k = 0.1210(0.0010)0.1240. One of the purpose of our 
calculation is to examine how the results depend on the 
choice of interpolating field operators. We employ several 
types of interpolating fields as (a) the NK*-type, (b) the 
(color-) twisted NK*-type, (c) a diquark-type, and adopt 
a smeared source to enhance the low-lying spectra. 

In J p = 3/2~ channel, we obtain massive states 
7715Q ~ 2.1 — 2.2 GeV except for the diquark-type inter- 
polating field, which involves a considerable size of the 
statistical error. In J p = 3/2+ channel, we obtain more 
massive states 7715Q ~ 2.4 — 2.6 GeV. None of these 5Q 
states appear below the NK threshold. Note that the NK 
threshold is raised up by about 200 — 250 MeV due to 
the finite extent of the spatial lattice as L ~ 2.15 fm, 
from which we expect the penta-quark signal to appear 
below the (raised) NK threshold considering the empir- 
ical mass difference between N+K(1440) and 8 + (1540). 
To clarify whether our 5Q states are compact resonance 
or not, we perform an analysis with the hybrid boundary 
condition(HBC), which was recently proposed in Ref . |6ll| . 
HBC analysis indicates that no compact 5Q resonance is 
contained in our 5Q states both in J p = 3/2 ± channels. 

The paper is organized as follows. In Sect. |H] we dis- 
cuss the general formalisms. We begin by introducing 
several types of interpolating fields, determining their 
parity transformation properties. We next consider the 
temporal correlator and its spectral decomposition. We 



finally discuss the two-particle scattering states involved 
in 5Q spectra, and introduce the hybrid boundary con- 
dition (HBC) to examine whether a state of our concern 
is a compact resonance state or not. Sect. Illll is devoted 
to the brief descriptions of our lattice action and param- 
eters. In Sect. IIVI we present our numerical results for 
J p = 3/2 ± channels in the standard periodic boundary 
condition(PBC). We show 5Q correlators of various in- 
terpolating fields, i.e., the NK*-type, the (color-)twisted 
NK*-type, the diquark-type. In Sect. E] we attempt to 
determine whether these 5Q states are compact 5Q res- 
onance states or two-particle scattering states by using 
the HBC. In Sect. IVll we summarize our results. 

II. GENERAL FORMALISMS 

A. Interpolating fields 

We consider an iso-scalar in terp ola ting field of NK*- 
type in Rarita-Schwinger form |73t l74l l75l | as 

IpH = £ abc (u^Cj 5 d b ) U c ■ (Sdj^dd) (1) 

where /i denotes the Lorentz index, a—d refer to the color 
indices, and C — 7472 denotes the charge conjugation 
matrix. Unless otherwise indicated, the gamma matrices 
are represented in the Euclidean form given in Ref . |76p . 

We are also interested in the ( color- )twisted NK* -type 
interpolating field as 

ipn = tabc (u^C^d b ) u d ■ {sd1^d c ) (2) 

-Cafcc (ul Cj 5 d b ) d d ■ (s d 7 M U c ) , 

which is an extension to the one originally proposed in 
Ref.[H3| to study J p = 1/2 P 5Q states. It has a slightly 
more elaborate color-structure than Eq. suggesting 
somewhat stronger coupling to a genuine 5Q state, if it 
exists, than simple NK* states. 

Another interpolating fields of our possible interests 
are diquark-type interpolating fields such as 

1/V = tabctdeftcfg (u^Cj 5 d b ) (u P C '757 M tf e ) C~f 5 S g (3) 

which is an extension to the one proposed in Refs. |3rj 
IHsj . The first factor corresponds to the scalar diquark 
(color 3, / = 0, J p = + ), which is expected to play 
important roles in hadron physics The second fac- 
tor corresponds to the vector diquark (color 3, / = 0, 
J =1"). Note that, although the axial- vector diquark 
(color 3, / = 1, J p = 1 + ) is considered to play a more 
important role than the vector diquark, it cannot replace 
the vector diquark due to its iso-vector nature. Unless 
otherwise indicated, we refer to Eq. © as the "diquark- 
type" interpolating field. We can also consider another 
interpolating field of diquark-type as 

IpH = tabctdeftcfg (u^Cd b ) {u^C^^de) Cs g , (4) 
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which consists of the pseudo-scalar diquark (color 3, 
I = 0, J p = 0~) and the vector diquark. However, 
actual lattice QCD calculation shows that its correlator 
is afflicted with quite a huge statistical error. A possible 
reason could be attributed to the fact that both of these 
diquark fields do not survive the non-relativistic limit. 
Hence, we do not consider this interpolating field in this 
paper. 

Under the spatial reflection of the quark fields as 

q(r,x) ^ j 4 q(T,-x), (5) 
all of these interpolating fields transform as 

1pi(T, x) -> -Ji1pi(T, -x), (6) 

for i = 1,2,3. 

B. 5Q correlators and parity projection 

We consider the Euclidean temporal correlator as 

G M „ (r) = (t, x)^ (0,0)), (7) 

X 

where J2s P r °j ec t s the total 5Q momentum to zero. 
Since the spin 3/2 contribution from the temporal com- 
ponent of Rarita-Schwinger spinor vanishes in the rest 
frame, we can restrict ourselves to the spatial parts, i.e., 
[if = 1, 2, 3. Now, Eq. is decomposed in the following 
way: 

Gij {r) = Pg^G^M+pg^G^Cr), (8) 

where i, j = 1, 2, 3 denote the spatial part of the Lorentz 
indices, p( 3 / 2 ) and P^ 1 / 2 ' denote the projection matrices 
onto the spin 3/2 and 1/2 subspaces defined as 



,(3/2) 



(l/3)7i7j, 



(9) 



Pg /2) = (1/3)77,. 
They satisfy the following relations as 



p(3/2) p (3/2) 
p(l/2)p(l/2) 



,(3/2) 

ik 

,(1/2) 



(10) 



,(1/2) 



p(l/2)p(3/2) 



p(3/2)p(l/2) = Q 



Here, summations over repeated indices are understood. 
G< 3 / 2 )(r) and G (1 / 2 >(t) in Eq. denote the spin 3/2 
and 1/2 contributions to G(r), respectively, which can 
be derived by operating p( 3 / 2 ) and pf 1 / 2 ) on G(r), re- 
spectively. (In our practical lattice QCD calculation, 
we construct the Rarita-Schwinger correlator Gy (r) for 
i = 1, 2, 3 and j = 3(fixed), and multiply P^ 3 / 2 ) from the 
left to obtain G< 3 / 2 )(t).) 



In the asymptotic region (0 <C r <C N t ), contamina- 
tions of excited states are suppressed. Considering the 
parity transformation property Eq. iJBJ, G^ 3 / 2 ^(r) and 
G^ x / 2 \t) are expressed in this region as 



G^' 2 \t) (11) 
= P+{|A 3/2 -| 2 e -™3 /2 - + |A 3/2+ | 2 e - (A, - r)m 3/ 2+ } 

- P-{\\ m+ \ 2 e- Tm ^ + +|A 3/2 -| 2 e - (JVt - r)m ^-} 

G (l/2) (r) 

= P+{|A 1/2 -| 2 e-™i/'- +|A 1/2+ | 2 e ^^- T W 2+ } 

- P_{|A 1/2+ | 2 e - Tm V2 + +| Al/2 _|2 e -W-r )m 1/2 _ j 5 



where P± = (1 ±7i)/2 denote the projection matrices 
onto the "upper" and "lower" Dirac subspaces, respec- 
tively. m 3 / 2 ± and mi«± denote the lowest- lying masses 
in J p = 3/2 ± and l/2 ± channels, respectively. A 3 / 2 ± and 
Ai/2± represent the couplings to the interpolating field 
Eq. with J p = 3/2 ± and l/2 ± states, respectively. In 
Eq. (|llfl . we adopt the anti-periodic boundary condition 
along the temporal direction. A brief derivation of Eq. (JSJ 
and Eq. Ijllfl is presented in Appendix. [S] The forward 
propagation is dominant in the region < r < Nt/2, 
while the backward propagation is dominant in the re- 
gion Nt/2 < r < Nt. To separate the negative (positive) 
parity contribution, we restrict ourselves to the region 
< t < Nt/2, and examine the "upper" ("lower") Dirac 
component. 



C. Scattering states involved in 5Q spectrum 

We consider the (two-particle) scattering states in- 
volved in 5Q spectrum. For J p = 3/2 ± iso-scalar penta- 
quarks, NK and NK* scattering states play an important 
role. (AK does not couple to the iso-scalar channel.) 
These states are expressed as 

\N(p,s)K(-p)), \N(p,s)K*(-p,i)), (12) 

where s and i denote the spin of the nucleon and K* , and 
p denotes the spatial momentum allowed for a particular 
choice of the spatial boundary condition adopted. For 
instance, if these hadrons are subject to the spatially pe- 
riodic boundary condition, their momenta are quantized 
as 



Pi 



2niTr/L, € 



(13) 



where L denotes the spatial extent of the lattice. In con- 
trast, if they are subject to the spatially anti-periodic 
boundary condition, their momenta are quantized as 



Pi 



{2n % + 1)tt/L, n t € 



(14) 
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K 




0.1210 


0.1220 


0.1230 


0.1240 


emp. 


NK* (s-wave) 


PBC 


2.996 


2.815 


2.633 


2.445 


1.830 


NK* (p-wave) 


PBC 


3.222 


3.052 


2.883 


2.710 


2.163 


NK(p/d-wave) 


PBC 


2.987 


2.806 


2.624 


2.438 


1.865 


NK*(s/p-wave) 


HBC 


3.167 


2.995 


2.823 


2.647 


2.084 


NK(p/d-wave) 


HBC 


2.924 


2.739 


2.553 


2.363 


1.770 



TABLE I: Numerical values of NK and NK* thresholds for 
each hopping parameter k in the physical unit GeV in the 
spatial lattice of the size L ~ 2.15 fm in PBC and HBC. The 
rightmost column labeled as "emp." shows the thresholds 
for the physical values of N, K, and K* as mjv — 0.94 GeV, 
m K ~ 0.5 GeV, m K * ^ 0.89 GeV. 

We first perform the parity projections. The positive and 
the negative parity states are obtained in the following 
way: 

\NK(±)) (15) 
= \N(p, s)K(-p)) T \N(-p, s)K{p)) 

\NK*(±)) (16) 
= \N(p,s)K*(-p,i)) T \N(-p,s)K*(p,i)). 

Assuming that the interactions between N and K and 
between N and K* are weak, their energies are approxi- 
mated as 

E NK ~ ^ m 2 N + p 2 + m 2 K + p 2 (17) 
E N . K ~ ^m 2 N +jP + ^/m 2 K , (18) 

respectively. The scattering states which couple to J p — 
3/2^ penta-quarks are obtained as spin-3/2 projections 
of Eq. JT3J) and Eq. JTJJ. The d-wave NK states and the 
s-wave NK* states can couple to the J p = 3/2" channel, 
while the p-wave NK states and NK* states can couple 
to the J p = 3/2+ channel. 

The scattering states with vanishing spatial momen- 
tum p — are exceptional in the following sense. On the 
one hand, the positive parity states vanish, because the 
first terms coincides with the second terms in Eq. i|15|) 
and Eq. l|16fl in the right hand side. On the other hand, 
the negative parity states are constructed only from the 
spin degrees of freedom, i.e., the spin degrees of freedom 
of the nucleon in Eq. I|15fl . and the spin degrees of free- 
doms of the nucleon and K* in Eq. I|16l) . By counting the 
degeneracy of the resulting states, it is straightforward 
to see that no d-wave states are contained, i.e., Eq. i|15|) 
gives only s-wave NK states in J p — 1/2" channel, and 
that Eq. (|Td|) gives only s-wave NK* states in J p = 1/2" 
and 3/2~ channels. 

D. Hybrid boundary condition(HBC) 

In order to determine whether a state of our concern 
is a compact 5Q resonance state or a scattering state 



of two particles, we use two distinct spatial boundary 
conditions(BC), i.e., the (standard) periodic BC(PBC) 
and the hybrid BC(HBC), which is recently proposed in 
Ref.[6l|. In PBC, one imposes the spatially periodic BC 
on u,d and s-quarks. As a result, all the hadrons are 
subject to the periodic BC. In this case, due to Eq. i|13[l . 
all hadrons can take zero-momentum, and the smallest 
non-vanishing momentum p min is of the form as 

(±271-/^,0,0), (0,±2tt/L,0), (0,0,±27r/£), (19) 

which gives 

\p™ c \=2n/L. (20) 

On the other hand, in HBC, we impose the spatially anti- 
periodic BC on u and d-quarks, whereas the spatially 
periodic BC is imposed on s-quark. Since N(uud,udd), 
K(us, ds) and K*(us, ds) contain odd numbers of u and d 
quarks, they are subject to the anti-periodic BC. There- 
fore, due to Eq. (|14fl . N, K and K* cannot have a vanish- 
ing momentum in HBC. The smallest possible momen- 
tum p m i n is of the form as 

(±tt/L,±7t/L,±tt/L). (21) 

Hence, its norm |p m in| is expressed as 

W£? C \=V37r/L. (22) 

In contrast, Q + (uudds) is subject to the spatially peri- 
odic BC, since it contains even number of u and d quarks. 
Therefore, + can have the vanishing momentum. 

Switching from PBC, HBC affects the low- lying two- 
particle scattering spectrum. A drastic change is ex- 
pected in the s-wave NK* channel. In PBC, the energy 
of the lowest NK* state is given as 

E^iNK* (s-wave)) ~ m N + m K ,. (23) 

In contrast, in HBC, since both N and K* are required to 
have non- vanishing momenta |p m i n | = vStt/L, the energy 
of the lowest NK* state is raised up as 

E™ C ( NK* (s-wave)) (24) 

~ yJm 2 N + 3tt 2 /L 2 + ^m\, + 3tt 2 /L 2 . 

Note that the shift amounts typically to a few hundred 
McV for L ~ 2 fm. 

HBC affects NK(d-wave), NK(p-wave), NK* (p-wave) 
as well. However, these changes are not as drastic as 
that in NK* (s-wave), because they are induced by the 
minor change in the minimum momentum from |p m i n | = 
2tt/L to yfl-n/L. In PBC, the energies of the lowest two- 
particles states are expressed as 

£ m P n c ( A^(p/d-wave )) (25) 

~ \Jm 2 N + 4tt 2 /£ 2 + ^Jm 2 K + 4tt 2 /L 2 

£ m P n c ( A^*(p-wave) ) 
~ ym 2 , + 47r 2 /£ 2 + \Jm\, +4tt 2 /L 2 . 
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In HBC, they are shifted down as 



E, 



HBC 



(NK (p/d-wave)) 



(26) 



3tt 2 /L 2 + , m 2 K + 3w 2 /L 2 



E, 



HBC 



(iYiT(p-wave)) 



N 



■3n 2 /L 2 



i 2 K ,+3n 2 /L 2 . 



Numerical values of NK and NK* thresholds for each 
hopping parameter in spatial lattice of the size L ~ 2.15 
fm for both PBC and HBC are summarized in Table |U 

In contrast to the scattering states, HBC is not ex- 
pected to affect a compact 5Q resonance + so much. 
Since Q + (uudds) can have vanishing momentum also in 
HBC, the shift of the penta-quark mass tosq originates 
only from the change in its intrinsic structure. In this 
case, the shift is expected to be less significant than the 
shift induced by the kinematic reason as is the case in N, 
K, and K*. Now our way to find a compact 5Q resonance 
state is to seek for such a state which is not affected by 
HBC. 



III. LATTICE ACTIONS AND PARAMETERS 

To generate gauge field configurations, we use the stan- 
dard plaquette action on the anisotropic lattice of the size 
12 3 x 96 as 



J2 ReTr{l-P^(x)} (27) 



JL1. 

N C JG 

^- 7G ReTr{l-P i4 (z)}, 



cc,i<j<3 



x,i<3 



where P^(x) G SU(3) denotes the plaquette operator 
in the /i-i/-plane. The lattice parameter and the bare 
anisotropy parameter are fixed as [5 = 2N c /g 2 — 5.75 
and 7g = 3.2552, respectively. These values are de- 
termined to reproduce the renormalized anisotropy as 
£ = a s /at = 4 69]. Adopting the pseudo-heat-bath al- 
gorithm, we pick up gauge field configurations every 500 
sweeps after skipping 10,000 sweeps for the thcrmaliza- 
tion. We use totally 1000 gauge field configurations to 
construct the temporal correlators. Note that the high 
statistics of iV con f = 1000 is quite essential for our study, 
because the 5Q correlators for spin 3/2 states are found 
to be rather noisy. In fact, a preliminary analysis with 
less statistics -/Vc 0n f ~ 500 leads to a spurious resonance- 
like state 78] . The lattice spacing is determined from 
the static quark potential adopting the Sommer param- 



eter r-g = 395 MeV (r ~ 0.5 fm) as a s 



1.100(6) 



GeV (<z s ~ 0.18 fm). Note that the lattice size 12 3 x 96 
amounts to (2.15fm) 3 x (4.30fm) in the physical unit. 
We adopt the <3(a)-improved Wilson (clover) action on 



the anisotropic lattice for quark fields ip an d ip as [7(J 

S F = ^(x)K{x,y)fl>(y), (28) 



K 



(x, y) = 5 x . y - re t { (1 - 74 ) Ua(x) 5 x+ly 

+(1 + 74)^(^-4) S x _ ly } 

- re s ^{(r- 7i ) Ui{x) S x+ly 

+ {r + ll ) U }{x-l)5 x _ ly ) 

i i<j 

where k s and K t denote the hopping parameters for the 
spatial and the temporal directions, respectively. The 
field strength F^ v is defined through the standard clover- 
leaf-type construction, r denotes the Wilson parameter. 
ce and cb denote the clover coefficients. To achieve the 
tadpole improvement, the link variables are rescaled as 
Ui(x) — > Ui{x)/u s and Un(x) — > U4(x)/ut, where u s and 
Ut denote the mean-field values of the spatial and tem- 
poral link variables, respectively |z3,[^]. This is equiv- 
alent to the redefinition of the hopping parameters as 
the tadpole-improved ones (with tilde), i.e., k s = k s /u s 
and K t = k t /ut- The anisotropy parameter is defined 
as 7f = kt/ks, which coincides with the renormalized 
anisotropy £ = a s /at for sufficiently small quark mass at 
the tadpole-improved level |70|. For given k s , the four 
parameters r, ce, cb and K s /K t should be, in principle, 
tuned so that "Lorentz symmetry" holds up to discretiza- 
tion errors of 0(a 2 ). Here, r, ce and cb are fixed by 
adopting the tadpole improved tree-level values as 



1 



CE 



u s u 2 : 



1 

cb = — r. 



(29) 



Only the value of k s / 'n t (= jp ■ (us/ut)) is tuned nonper- 
turbatively by using the meson dispersion relation [70|. 
It is convenient to define k as 



11, 

- = ^- - 2 (7jr - 3r - 4) . 



(30) 



Then the bare quark mass is expressed as mo = — 
8) in the spatial lattice unit in the continuum limit. This 
k plays the role of the hopping parameter "re" in the 
isotropic formulation. For detail, see Refs. [Ulnl, where 
we take the lattice parameters. The values of the lattice 
parameters are summarized in Table llll 

We adopt four values of the hopping parameter as 
re = 0.1210,0.1220,0.1230 and 0.1240, which correspond 
to m^/nip = 0.81, 0.78, 0.73 and 0.66, respectively. These 
values roughly cover the region m s < m < 2m s . For tem- 
poral direction, we impose anti-periodic boundary condi- 
tion on all the quark fields. For spatial directions, we 
impose periodic boundary condition on all the quarks, 
unless otherwise indicated. We refer to this boundary 
condition as "periodic BC (PBC)". 



6 



P 7G a s /at a a 1 [GeV] Size N con i u s ut 7f k c Values of k 

5.75 3.2552 4 1.100(6) 12 3 X 96 1000 0.7620(2) 0.9871(0) 3.909 0.12640(5) 0.1240, 0.1230, 0.1220, 0.1210 



TABLE II: Parameters of the lattice simulation. The spatial lattice spacing a a is determined with r^ 1 — 385 MeV for the 
Sommer parameter. The mean-field values of link variables (w CT and u T ) are defined in the Landau gauge. n c denotes the critical 
value of «. 



0.1210 



0.1220 



0.1230 



0.1240 



^phys. 



m p 
m K 
m K * 



1.007(2) 
1.240(4) 
0.846(2) 
1.119(6) 
1.877(4) 



0.897(1) 
1.157(5) 
0.785(1) 
1.076(7) 
1.739(3) 



0.785(2) 
1.074(7) 
0.722(2) 
1.033(9) 
1.600(4) 



m N * 2.325(17) 2.194(21) 2.059(28) 



0.658(2) 

0.991(11) 

0.658(2) 

0.991(11) 

1.454(5) 

1.918(42) 



0.140 

0.823(13) 

0.476(2) 

0.902(15) 

1.164(8) 

1.648(53) 



TABLE III: Masses of tt, p, K, K*, N, and N* for each hopping 
parameter k in the physical unit GeV. « p hys. — 0.1261 denotes 
the value of k which achieves m v ~ 0.14 GeV. 



By keeping k s = 0.1240 fixed for s quark, and by 
changing k — 0.1210 — 0.1240 for u and d quarks, we 
perform the chiral extrapolation to the physical quark 
mass region. In the following part of the paper, we will 
use 



(k s ,k) = (0.1240,0.1220), 



(31) 



as a typical set of hopping parameters in presenting cor- 
relators and effective mass plots. For convenience, we 
summarize masses of tt, p, K, K*, N and N*(J P = l/2~ 
baryon) for each hopping parameter k together with their 
values at the physical quark mass in Table HTTI Here, the 
chiral extrapolations of these particles are performed by 
a linear function in m 2 .. Unless otherwise indicated, we 
adopt the jackknife prescription to estimate the statisti- 
cal errors. 

We use a smeared source to enhance the low-lying spec- 
tra. More precisely, we employ spatially extended inter- 
polating fields of the gaussian size p ~ 0.4 fm, which is 
obtained by replacing the quark fields q(x) in 5Q inter- 
polating fields by the smeared quark fields q S mcar(a;) in 
the Coulomb gauge as 



r (r,f) = AA^cxp^ - 



P - y\ 

2p 2 



Q(r,ff), (32) 



where Af is an appropriate normalization factor. For a 
practical use, we extend Eq. Ij32(l appropriately so as to fit 
a particular choice of the spatial boundary condition. In 
this paper, we present correlators with a smeared source 
and a point sink. 



IV. NUMERICAL RESULTS ON 5Q SPECTRUM 

We present our lattice QCD results on 5Q spectrum in 
the standard periodic boundary condition(PBC) in this 



section. 



A. J p = 3/2' 5Q spectrum in PBC 

We consider 5Q spectrum in J p = 3/2~ channel. In 
Fig. ^ we show the effective mass plots in J p = 3/2 ~~ 
channel for three interpolating fields, i.e., (a) the NK*- 
type, (b) the twisted NK*-type, (c) a diquark-type. The 
dotted lines indicate the s-wave NK* and the d-wave 
NK thresholds, which happen to coincide accidentally in 
Fig. fflfor the spatial lattice size L ~ 2.15 fm. 

We define the effective mass as a function of r by 

™Mr) = log ( ^(3/ 2 ) (r+r J > ( 33 ) 

where G^ 3 ^ 2 ^(t) denotes the temporal correlator. At 
sufficiently large r, the correlator is dominated by the 
lowest-lying state with energy m as G^ 3 ' 2 '(t) ~ Ae~ mT . 
Then Eq. (lo'.'A) gives a constant as m e ff(r) ~ m. Thus a 
plateau in m e ff(r) indicates that the correlator is satu- 
rated by a single-state. In such cases, we can perform a 
single-exponential fit in the plateau region. 

Fig. ^ ( a ) shows the effective mass plot for the NIC- 
type interpolating field. In the region < r < 24, the 
contamination of the higher spectral contributions are 
gradually reduced, which is indicated by the decreases in 
Tn e s(r). There is a plateau in the interval 25 < t < 35, 
where a single-state is expected to dominate the 5Q corre- 
lator. Beyond r ~ 36, the statistical error becomes large. 
In addition, the effect of the backward propagation be- 
comes gradually more significant as r ~ 48 is approached. 
Hence, we simply neglect the data for r > 36, and per- 
form the single-exponential fit in the region 25 < r < 35. 
We obtain m5Q = 2.90(2) GeV, which is denoted by the 
solid line. One sees that the 5Q states appears above the 
s-wave NK* and the d-wave NK thresholds. 

Fig.n(b) shows the effective mass plot for the twisted 
NK*-type interpolating field. There is a plateau in the 
interval 24 < r < 35, where the single-exponential fit 
is performed leading to m 5 q = 2.89(1) GeV. The 5Q 
state is again above the s-wave NK* and the d-wave NK 
thresholds. 

Fig. [I](c) shows the effective mass plot for the diquark- 
type interpolating field. We see that the statistical er- 
ror is too large to identify the plateau unambiguously. 
Hence, we do not perform the fit. Note that this plot is 
obtained by using iV con f = 1000 gauge configurations. A 
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(a) NK* (J =3/2 ) [PBC] 



J =3 /2~ 



NK* (s-wave) 
NK (d-wave) 
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t [a t ] 

(b) twisted NK* (J P =3/2~) [PBC] 



NK* (s-wave) 
NK (d-wave ) 
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t [a t ] 

(c) diquark ( J P =3/2~) [PBC] 
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FIG. 1: The 5Q effective mass plots in J p — 3/2 ~ channel in 
the standard periodic BC(PBC) for three types of interpolat- 
ing fields, i.e., (a) the NK*-type, (b) the twisted NK*-type, 
and (c) the diquark-type. Eq. 13111 is adopted as a typical 
set of hopping parameters. The statistical error is estimated 
with the jackknife prescription. The dotted lines indicate the 
s-wave NK* and the d-wave NK threshold in the spatial lat- 
tice size L ~ 2.15 fm. Note that they accidentally coincide 
with each other. The solid lines denote the results of the 
single-exponential fit performed in each plateau region. 



3.5 



> 

QJ 
U 
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2.5 
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[GeV 



FIG. 2: msQ in J p = 3/2 _ channel against m\ for the two in- 
terpolating fields, i.e., (circle) the NK*-type, (box) the twisted 
NK*-type. Open symbols denote the direct lattice QCD data, 
while the closed symbols and the solid lines represent the re- 
sults of the chiral extrapolations to the physical quark mass 
region. 



J p 


I.F. 


0.1210 


0.1220 


0.1230 


0.1240 


^phys. 


3/2- 


(a) 


3.08(1) 


2.90(2) 


2.72(2) 


2.54(3) 


2.17(4) 


3/2" 


(b) 


3.08(1) 


2.89(1) 


2.70(2) 


2.49(3) 


2.11(4) 


3/2" 


(c) 












3/2+ 


(a) 


3.52(2) 


3.34(3) 


3.17(11) 


3.00(5) 


2.64(7) 


3/2+ 


(b) 


3.27(3) 


3.11(4) 


2.95(5) 


2.83(9) 


2.48(10) 


3/2+ 


(c) 


3.34(2) 


3.16(2) 


2.98(3) 


2.78(5) 


2.42(6) 



TABLE IV: msQ for each value of n in the physical unit GeV. 
The first column labeled by "J p " indicates the spin and the 
parity. The second column labeled by "I.F." indicates the 
interpolating field used, i.e., (a) the NK*-type, (b) the twisted 
NK*-type, (c) the diquark-type. «;ph ys . — 0.1261 denotes the 
value of k which achieve ~ 0.14 GeV. "-" indicates that 
the fitting is not performed due to the large statistical error. 



possible reason for such a large noise is that the interpo- 
lating field Eq. © does not survive the non-relativistic 
limit due to the vector diquark. 

Now, we perform the chiral extrapolation. As men- 
tioned before, we keep k = 0.1240 fixed for s-quark, and 
vary k = 0.1210 — 0.1240 for u and d quarks. Fig.|2|shows 
the 5Q masses in J p = 3/2 _ channel against m^. Circles 
and boxes denote the data obtained from the NK*-type 
and the twisted NK*-type 5Q correlators, respectively. 
Note that they agree with each other within the statis- 
tical error. The open symbols refer to the direct lattice 
QCD data. Since these data behave almost linearly in 
mj, we adopt the linear chiral extrapolation in to 
obtain tosq in the physical quark mass region. Note that 
the ordinary non-PS mesons and baryons show similar 
linearity in [7l|. The closed symbols denote the re- 
sults of the chiral extrapolation. We see that all the 5Q 
states appear above the s-wave NK* and the d-wave NK 
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thresholds. As a result of the chiral extrapolation, we 
obtain only massive 5Q states as m 5 Q = 2.17(4), 2.11(4) 
GeV from the NK*-type and the twisted NK*-type cor- 
relators, respectively, which is too heavy to be identified 
with the experimentally observed + (154O). Numerical 
values of ttisq at each hopping parameter together with 
their chirally extrapolated values are summarized in Ta- 
ble ]W\ To obtain a low- lying state at msQ ~ 1540 MeV, 
a 5Q state should appear below these thresholds at least 
in the light quark mass region. In this case, a significantly 
large chiral effect is required. Of course, this point can 
be in principle clarified by an explicit lattice QCD calcu- 
lation with chiral fcrmions. 



B. J p = 3/2+ 5Q spectrum in PBC 

We consider 5Q spectrum in J p = 3/2 + channel. 
J p = 3/2 + is an interesting quantum number from the 
view point of the diquark picture of Jaffe and Wilczek[22] • 
In this picture, the pair of diquarks has angular mo- 
mentum one, which is combined with the spin 1/2 of 
s quark. Hence, there are two possibilities as J p = 
l/2 + and 3/2 + , i.e., the diquark picture can support 
J p = 3/2 + possibility as well. Its mass splits from the 
J p = 1/2+ state depending on a particular form of the 
LS-interaction. If it is massive, it is expected to have a 
large decay width. If it is light enough, its exotic struc- 
ture may work to implement the narrow decay width as 
in J p = 1/2+ case. 

In Fig. |3 we show the 5Q effective mass plots in PBC 
employing three types of interpolating fields, i.e., (a) the 
NK*-type, (b) the twisted NK*-type, (c) the diquark- 
type. The dotted lines indicate the s-wave N*K*, the 
p-wave NK* and the p-wave NK thresholds in the spatial 
lattice of the size L ~ 2.15 fm, respectively. 

Fig. |2| (a) shows the 5Q effective mass plot employ- 
ing the NK*-type interpolating field. In the region, 
< t < 17, the contaminations of higher spectral con- 
tributions become gradually reduced. There is a flat re- 
gion 18 < t < 30, which is still afflicted with slightly 
large statistical errors. The single-exponential fit in this 
region gives m 5 Q = 3.34(3) GeV. Note that this value 
agrees with the s-wave N*K* threshold E th ~ 3.27 GeV. 
(See Table Hill for m N ,.) 

Fig. |3|(b) shows the 5Q effective mass plot correspond- 
ing to the twisted NK*-type interpolating field. We have 
a rather stable plateau in the interval 21 < r < 27, 
where the single-exponential fit is performed. We ob- 
tain 7715Q = 3.11(4) GeV. The result is denoted by the 
solid line. 

Fig. 13 (c) shows the 5Q effective mass plot for the 
diquark-type interpolating field. We find a plateau in 
the interval 19 < r < 29, where the single-exponential fit 
is performed. We obtain m^Q = 3.16(2) GeV, which is 
denoted by the solid line. 

Now, we perform the chiral extrapolation. In Fig. 0] 
TO 5Q is plotted against m„. Circles, boxes and trian- 
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(b) twisted NK* (J P =3/2 + ) [PBC] 
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(c) diquark (J P =3/2 + ) [PBC] 
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FIG. 3: The 5Q effective mass plots in J p — 3/2 + channel 
in PBC for three types of interpolating fields, i.e., (a) the 
NK*-type, (b) the twisted NK*-type, and (c) the diquark- 
type. Eq. I31H is adopted as a typical set of hopping parame- 
ters. The dotted lines indicate the s-wave N*K*, the p-wave 
NK* and the p-wave NK thresholds in the spatial lattice size 
L ~ 2.15 fm. The solid lines denote the results of the single- 
exponential fit performed in each plateau region. 
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V. ANALYSIS WITH HBC 

In the previous section, we have only massive 5Q 
states, which are obtained by using the linear chiral ex- 
trapolation in mj. However, the chiral behavior may 
deviate from a simple linear one in the light quark mass 
region, which could lead to somewhat less massive states. 
Considering this, we think it of worth at this stage to an- 
alyze whether our 5Q states are compact 5Q resonances 
or not. This is done by switching the spatial periodic BC 
to the hybrid BC(HBC) introduced in Sect. OH 



FIG. 4: ni5Q in J p — 3/2 + channel against m\ for the 
three interpolating fields, i.e., (circle) the NK*-type, (box) 
the twisted NK*-type, and (triangle) a diquark-type [Eq. J3J]. 
Open symbols denote the direct lattice QCD data, while the 
closed symbols and the solid lines represent the results of the 
chiral extrapolations to the physical quark mass region. 



gles denote the data obtained from the NK*-type, the 
twisted NK*-type and the diquark-type 5Q correlators, 
respectively. Note that the latter two agree with each 
other within the statistical errors. 

As a result of the chiral extrapolation, we obtain 
m 5Q = 2.64(7) GeV from the NK*-type correlator, 
to 5 q = 2.48(10) GeV from the twisted NK*-type corre- 
lator, and m 5 Q = 2.42(6) GeV from the diquark-type 
correlator. Numerical values of m$Q in J p = 3/2+ 
channel at each hopping parameter together with their 
chirally extrapolated values are also summarized in Ta- 
bic IIVI The two data from the twisted NK*-type and 
the diquark-type correlators are considered to be almost 
consistent with the p-wave NK* threshold, while the data 
from the NK*-type correlator seems to correspond to 
a more massive state, which is likely to be consistent 
with the N*K* (s-wave) threshold. We see again that all 
of our data of m^Q appear above the NK* threshold(p- 
wave), which is located above the artificially raised NK 
threshold(p-wave) due to the finiteness of the spatial lat- 
tice as L ~ 2.15 fm. As a result, we are left only with 
such massive 5Q states. 

Now, several comments are in order. (1) Ref.[62j re- 
ported the existence of a low-lying 5Q state in J p = 3/2+ 
channel using NK*-typc interpolating field. However, we 
have not observed such a low-lying 5Q state in our calcu- 
lation. (2) Recall that, except for a single calculation |59|. 
lattice QCD calculations indicate that J p = l/2 + state 
is heavy [57l l58l l60l l6ll I62L lo3L l64l l65l l66| . for instance 
to 5 q ~ 2.25 GeV in Ref. 61]. From the viewpoint of the 
diquark picture, it could be natural to obtain such mas- 
sive 5Q states in J p = 3/2+ channel. If there were a 
low- lying 5Q state in J p = 3/2 + channel, then the di- 
quark picture could suggest also a low-lying 5Q state in 



A. J p = 3/2 5Q spectrum in HBC 

Fig. [5] shows the 5Q effective mass plots in HBC em- 
ploying the three types of interpolating fields, i.e., (a) 
the NK*-type, (b) the twisted NK*-type, and (c) the di- 
quark type. These figures should be compared with their 
PBC counterparts in Fig. ^ The dotted lines denote 
the s-wave NK* and the d-wave NK thresholds. For the 
typical set of hopping parameters, i.e., Eq. the s- 

wave NK* threshold(the thick dotted line) is raised up 
by ~ 180 MeV, and the d-wave NK threshold(the thin 
dotted line) is lowered down by ~ 70 MeV due to HBC 
in the finite spatial extent as L ~ 2.15 fm. (See Tabled) 

Fig. (a) shows the 5Q effective mass plot for the 
NK*-type interpolating field in HBC. We find a plateau 
in the interval 23 < r < 35, where the single-exponential 
fit is performed leading to 7715Q — 2.98(1) GeV, which is 
denoted by the solid line. We see that 7715Q is raised up 
by 80 MeV due to HBC. The value of msQ is consistent 
with the s-wave NK* threshold within the statistical er- 
ror. Therefore, we regard this state as an NK* scattering 
state. 

Fig- HI (b) shows the 5Q effective mass plot for the 
twisted NK*-type interpolating field. We find a plateau 
in the interval 24 < r < 35, where the single exponential 
fit is performed leading to m^Q = 2.98(1) GeV, which is 
denoted by the solid line. The situation is similar to the 
NK*-interpolating field case. We see that m^Q is raised 
up by 90 MeV due to HBC. Since the value is consis- 
tent with the s-wave NK* threshold within the statistical 
error, we regard it as an NK* scattering state. 

Fig. (c) shows the 5Q effective mass plot for the 
diquark-type interpolating field. We see that it is afflicted 
with considerable size of statistical errors as before, due 
to which the best-fit is not performed. 

In this way, all of our 5Q states in J p = 3/2~ channel 
turn out to be NK* scattering states. More precisely, we 
do not observe any compact 5Q resonance states in J p = 
3/2~ channel below the raised s-wave NK* threshold, i.e., 
in the following region: 

E < ^m 2 N +pL in + \jm\, +^ min , (34) 
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(c) diquark (J P =3/2~) [HBC] 
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FIG. 5: The 5Q effective mass plots in J p = 3/2 ~ channel in 
HBC for three types of interpolating fields, i.e., (a) the NK*- 
type, (b) the twisted NK*-type, and (c) the diquark-type. 
The meanings of the dotted lines and the solid lines are the 
same as in Fig. 
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(b) twisted NK* (J P =3/2 + ) [HBC] 
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(c) diquark ( J P =3/2 + ) [HBC] 
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FIG. 6: 5Q effective mass plots in J p = 3/2+ channel in HBC 
for three types of interpolating fields (a) the NK*-type, (b) the 
twisted NK*-type, and (c) the diquark-type. The meanings 
of the dotted and the solid lines are the same as in Fig. [3] 



Fig. shows the 5Q effective mass plots in HBC em- 
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ploying the three types of interpolating fields, i.e., (a) 
the NK*-type, (b) the twisted NK*-type, and (c) the 
diquark-type. These figures should be compared with 
their PBC counterparts in Fig. [21 The meanings of the 
dotted and the solid lines are the same as in Fig. [21 

HBC may not be useful in J p — 3/2 + channel, since it 
induces only minor changes in the two-particle spectra. 
For the typical set of hopping parameters, i.e., Eq. (|31|l . 
the p-wave NK* threshold is lowered down only by ~ 
60 MeV, and the p-wave NK threshold is lowered down 
only by ~ 70 MeV. We see that these shifts are rather 
small. This is because they are induced by the changes 
in the minimum non-vanishing momentum, i.e., |p m i n | = 
2tt/L ~ 576 MeV to y/Zn/L ~ 499 MeV as mentioned 
before. In J p = 3/2+ channel, N*K*(s-wave) threshold 
shows the most drastic change, i.e., the upper shift by 170 
MeV, which however plays a less significant role, since its 
location is at rather high energy. 

Fig. [HI (a) shows the 5Q effective mass plot employing 
the NK*-type interpolating field. There is a flat region 
16 < t < 25, which is still afflicted with slightly large 
statistical errors. The single-exponential fit in this region 
leads to m 5 Q = 3.38(2) GeV, which is denoted by the 
solid line. We see that m 5 Q is raised up by 40 MeV. 
Although the shift of 40 MeV is rather small, m^Q is 
again almost consistent with the s-wave N*K* threshold. 
Considering its rather large statistical error, this 5Q state 
is likely to be an s-wave N*K* scattering state. To draw 
a more solid conclusion on this state, it is necessary to 
improve the statistics further more. 

Fig. El (b) shows the 5Q effective mass plot employ- 
ing the twisted NK*-type interpolating field. There is a 
plateau in the interval 23 < r < 31, where we perform 
the single-exponential fit. The result m^q — 3.02(3) GeV 
is denoted by the solid line. We see that m^Q is lowered 
down by 90 MeV, which is considered to be consistent 
with the shift of the NK* (p-wave) threshold. Therefore, 
this state is likely to be an NK* (p-wave) scattering state. 

Fig. (c) shows the 5Q effective mass plot employing 
the diquark-type interpolating field. Although the data 
is slightly noisy, there is a plateau in the interval 23 < 
t < 28. A single-exponential fit in this plateau region 
leads to m^Q — 3.08(4) GeV, which is denoted by the 
solid line. m^Q is lowered down by 80 MeV due to HBC. 
The situation is similar to Fig. [51(b). This state is likely 
to be an NK* (p-wave) scattering state. 

In this way, all of our 5Q states are likely to be ei- 
ther N*K* (s-wave) or NK* (p-wave) states rather than 
compact 5Q resonance states. Of course, because HBC 
induces only minor changes in the 5Q spectrum in J p = 
3/2+ channel, and also because 5Q correlators still in- 
volve considerable size of statistical error, more statistics 
is desirable to draw a more solid conclusion on the real 
nature of these 5Q states. Here, we can at least state 
that these 5Q states are all massive, which locate above 
NK* (p-wave) threshold. 



VI. SUMMARY AND CONCLUSION 

We have studied J p = 3/2 ± penta-quark(5Q) baryons 
in anisotropic lattice QCD at the quenched level with 
a large number of gauge field configurations as N con { — 
1000 for high precision measurements. We emphasize 
that the spin of O + (1540) has not yet determined exper- 
imentally, and that the J p = 3/2~ assignment provides 
us with one of the possible solutions to the puzzle of the 
narrow decay width of 9 + (1540) [4^. We have employed 
the standard Wilson gauge action on the anisotropic lat- 
tice of the size 12 3 x 96 with the renormalized anisotropy 
fls/flt = 4 at j3 = 5.75, which leads to a s ~ 0.18 fm 
and at — 0.045 fm. We have found that correlators of 
5Q baryons in J p = 3/2 ± channels are rather noisy. 
Hence, the large statistics as iV con f = 1000 has played 
a key role to get a solid result in our calculation. For 
the quark part, we have employed 0(a)-improved Wilson 
(clover) action with four values of the hopping parame- 
ters as k — 0.1210(0.0010)0.1240, which roughly cover 
the quark mass region as m s < m < 2m s . To avoid 
the contaminations of higher spectral contributions, we 
have employed the spatially extended source in the 5Q 
correlators. 

We have examined several types of the 5Q interpolat- 
ing fields as (a) the NK*-type, (b) the (color-)twisted 
NK*-type, (c) the diquark-type. In J p = 3/2~ chan- 
nel, there are plateaus in the effective mass plots for the 
NK*-type and the twisted NK*-type interpolating field, 
whereas no plateau has been identified in that for the 
diquark-type interpolating field due to the significantly 
large statistical error. The former two give almost identi- 
cal results. We have employed the linear chiral extrapola- 
tions in m%, which have lead to m^Q ~ 2.17 and 2.11 GeV 
for the NK*-type and the twisted NK*-type 5Q correla- 
tors, respectively. In J p = 3/2+ channel, we have recog- 
nized plateaus in all the three effective mass plots. How- 
ever, the plateau for the NK*-type interpolating field is 
located at a somewhat higher energy than the other two. 
The chiral extrapolations have lead to m^Q ~ 2.64 GeV 
for the NK*-type correlator, m 5 Q ~ 2.48 GeV for the 
twisted NK*-type correlator, and m 5 Q ~ 2.42 GeV for 
the diquark-type correlator. In this way, our data have 
not supported low- lying 5Q states in both J p — 3/2 ± 
channels. All the 5Q states have been observed to appear 
above the d/p-wave NK threshold, which is artificially 
raised up by a few hundred MeV due to the finiteness of 
the spatial lattice as L ~ 2.15 fm. Note that, to obtain 
low- lying 5Q states in J p = 3/2 ± channel, a 5Q state 
should appear below the raised NK threshold (p/d- wave) 
at least in the light quark mass region. 

In order to clarify whether the observed states are com- 
pact 5Q resonances or not, we have performed an anal- 
ysis with hybrid boundary condition(HBC), which was 
recently proposed by Ref.[6ll|. In J p = 3/2~ channel, 
our 5Q states observed in the NK*-type and the twisted 
NK*-type correlators have turned out to be s-wave NK* 
scattering states. In J p = 3/2+ channel, for the twisted 
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NK*-type and the diquark-type correlators, the observed 
5Q states are most likely to be NK*(p-wave) scattering 
states. For the other one, i.e., the NK*-type interpolating 
field, although more statistics is needed to draw a definite 
conclusion, it is most likely to be an s-wave N*K* scat- 
tering state. Note that, since HBC does not affected the 
two-particle spectra so much in J p = 3/2+ channel, it is 
not easy to elucidate the natures of the 5Q states only 
with HBC. At any rate, whatever the real nature of these 
5Q states may be, they result in a considerably massive 
states in the physical quark mass region, which cannot be 
identified as 9 + (1540) without involving a significantly 
large chiral contribution. 

In this way, we have not obtained any relevant sig- 
nals for low-lying compact 5Q resonance states in J p = 
3/2^ channel below 2.1 GeV in this paper, although the 
J p = 3/2~ possibility provides us with one of the pos- 
sible solutions to the puzzle of the narrow decay width 
of 8 + (1540). To get more solid conclusion on the pen- 
taquark, it is important to perform the systematic stud- 
ies of the 5Q states with the various quantum numbers 
in lattice QCD with more sophisticated conditions. For 
instance, it is desired to use (1) unquenched full lat- 
tice QCD, (2) finer and larger volume lattice, (3) chiral 
fermion with small mass, (4) more sophisticated inter- 
polating field corresponding to the diquark picture and 
so on. In any case, the mysterious exotic hadron of the 
pentaquark would be much clarified in future studies of 
lattice QCD as well as in the future experiments. 
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APPENDIX A: SPECTRAL REPRESENTATION 

Considering the importance of the parity determina- 
tion of + , we present a brief derivation of the spec- 
tral representation of Rarita-Schwinger correlators, i.e., 
Eq. (JSJ, with Eq. Qll|) . In this section, gamma matri- 
ces are represented in Minkowskian form (See Ref.[7jJ). 
To avoid unnecessary complexities, we derive only J p = 
3/2± parts. The J p = 1/2* parts can be obtained as a 
slight modification, on which we will make a comment at 
the end of the section. 

We first consider the coupling of our interpolating 
fields ip^ to J p = 3/2 ± (anti-)baryon states. Due to 
Eq. JfjJ), our interpolating fields, i.e., Eqs. CQ,© and Q 
have the negative intrinsic parity. Hence, their couplings 



to J p = 3/2 (anti-)baryons are parameterized in the 
following way: 

(0\tp l i(0)\B 3 / 2 - (k, a)) = X 3/2 -u f _ l (m 3/ 2-;k,a)(Al) 
(01^(0)153/2- {k, a)) = Xl /2 -v f _ l (m 3/2 -;k,a), 

where \B 3 / 2 - (k, a)} and \B 3 / 2 - (k, a)} denote J p = 3/2~ 
(anti-)baryon states with momentum fc, helicity a, and 
mass m 3 / 2 +- u^(m]k,a) and u M (m;/c,o;) denote the 
Rarita-Schwinger spinors for J = 3/2 particles with 
momentum k, helicity a and mass m [TH. \l\ l75| . 
Eqs. CEJ,© and (2J couples to J p = 3/2+ (anti-)baryons 
as well. In this case, their couplings involve 75 in the fol- 
lowing way: 

(0|V^(0)|B 3/2 +(fc,a)) = X 3/2 +"f 5 u fl (m 3/2 +;k,a)(A2) 
(0|Vv(0)|S 3/2 +(fc,a;)) = - A 3/2+^( m 3/2+;^,a)75, 

where \B 3 / 2 + (k, a)) and \B 3 / 2 + {k, a)} denote J p — 3/2+ 
(anti-)baryon states with momentum fc, helicity a and 
mass TO3/2+ ■ " — " originates from the anti-commutativity 
of 70 and 75. 

To derive the spectral representation, the best way 
would be to express it in the operator representation in 
the following way: 

G^{t,x) = Z-'Tr (e- /?H T r [^(r,af)Vv(0)]) , (A3) 

where denotes the temporal extent of the lattice, H = 
Hq C d denotes the QCD Hamiltonian, Z = Tr(e _/3ff ) 
denotes the partition function, and T T [*] represents the 
time-ordered product along the imaginary time direc- 
tion. The interpolating fields are represented in the 
Heisenberg picture in imaginary-time, i.e., "0^ (t, 2;) = 
e Tff VM(0,x)e- rff and Vv(r,£) = e TH tP tl (0, x)e~ TH . By 
restricting ourselves to the interval < r < /3, Eq. (|A3|) 
reduces to 

G^(t,x) = Tr ( -_VyMVv(0) J . (A4) 

Note that it can be equivalently expressed as 

G^(r, x) = Tr Ufy ( T - 0, f)^W(0) J • (A5) 

In the large (3 limit, we have e~@ H jZ ~ |0)(0|, which 
is inserted into Eq. (|A4|) and Eq. I|A5|) . Note that the 
resulting two expressions serve as independent contribu- 
tions to the original "Tr", i.e., Eq. (jA4|) (or Eq. (|A5|) 1. 
Hence, we keep these two contributions to obtain 

G^{t,x) (A6) 
~ (0|Vy(T,f)Vv(0)|0) + (0|Vv(0)<V(t-/?,x)|0>. 

Note that the 1st term corresponds to the forward prop- 
agation, whereas the 2nd term to the backward prop- 
agation. By inserting single- (anti-)baryon intermediate 
states, and by using Eq. ijAlfl and Eq. (|A2|) . we are left 
with 
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4 r 

Ghv{t,x) = X/ y 



d 3 k m 3/2+ . 
(2nf k 



(0|?/y (f)|B 3/a+ (fc, a)) (B 3/2+ (k, a)|^(0)|0) 



(A7) 



3/2- 



3/2-* 



<i 3 fc m 3 / 2 - 

d 3 k m 3/2 + 
(2^) 3 fco 



(-l)^ /2) (fc) 



g-rfcogifc-x 


f m 3 / 2 - + 






V 2m 3/2" 






' "13/2+ - 






V 2^3/2+ 





i g-(/5-r)fc g-ife-S 



e -{p-T)k 0e -ik,x 



2m 



3/2- 



m 3/2 + +^ 



2m 



3/24 



where fco 



k 2 for J p 



3/2-, fc 



m- 



"3/2+ 

tities are used 



3/2- 

fc 2 for J p = 3/2+ , and the following iden- 



4 

V W ,(m;t,«KKi,«) = ~ ! ^^ /2) MA8) 

^ — ' 2m p 

Q = l 

4 m-fc 

^2v ll {m;k,a)v v {m;k,a) = — P^J 2) {k), 



( 3 /2) 

where P/V (fc) is the spin 3/2 projection operator de- 
fined as 

^ 3 / 2) (*0 = 9»u - - ^ #7/**" + *m7.>*0 • (A9) 

By performing the integration over x for zero-momentum 
projection, and by replacing the Minkowskian gamma 
matrices by their Euclidean counterparts, we finally ar- 
rive at the spectral representation (Eq. (JHJ with Eq. Ulljl). 

The derivation of the spin 1/2 parts is obtained by 
repeating a similar procedure using the following param- 



eterizations instead of Eq. ijAlf) and Eq. (|A2|) as 

(Ol^Wlfli/a- (fc, a)) (A10) 

= (Ai/ 2 -7m + A i/2- fc p) «(mi/ 2 -;fc,a) 
(0fe(0)|B 1/2 ~(fc,a)) 

= v(m 1/2 -;k,a) (a* /2 _7 jU - A'^-fc^J 
(01^(0)1^/2+ (fc, a)) 

= (Ai/ 2 +7p + Ai/2- fc/ij j5u(m 1/2 + ; fc, a) 
(0|^(0)|S 1/2 -(fc,a)) 

= -v{m x /2- ; fc, a)75 (a i/2 -7m - A'*/ 2 -fcJ , 

where |Bi/ 2 ±(fc, a)) and |Pi/2± (fc, a)) denote the J p = 
l/2 ± (anti-)baryon states with momentum fc, helicity a 
and mass mi/ 2 ±. u{m;k,a) and v(m;fc,a) denote the 
Dirac bispinors for spin 1/2 particles with mass m, mo- 
mentum fc and helicity a. \\/2± and A' 1 y 2± represent the 

couplings to J p = l/2 ± (anti-)baryons. 
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